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Abstract 

We study D-dimensional polymerized membranes embedded in d dimensions using a self- 
consistent screening approximation. It is exact for large d to order 1/d, for any d to order 
e = 4 — D and for d = D. For flat physical membranes (D — 2, d — 3) it predicts a roughness 
exponent ( = 0.590. For phantom membranes at the crumpling transition the size exponent is 
v = 0.732. It yields identical lower critical dimension for the flat phase and crumpling transition 
Dic(d) = Jn; (Di c = \[2 for codimension 1). For physical membranes with random quenched 
curvature ( = 0.775 in the new T = flat phase in good agreement with simulations. 

PACS: 64.60Fr,05.40,82.65Dp 



1 



There are now several experimental realizations of polymerized or solid-like membranes, 
such as protein networks of biological membranes 1 ' 2 , polymerized lipid bilayers 3 and some 
inorganic surfaces 4 . Unlike linear polymers, two dimensional sheets of molecules with fixed 
connectivity and non zero shear modulus are predicted to exhibit a flat phase with broken 
orientational symmetry Out of plane thermal undulations of solid membranes which induce 
a non-zero local Gaussian curvature are strongly suppressed because they are accompanied by 
in-plane shear deformations 5 . As a result, even "phantom" tethered membranes should be flat 
at low temperatures 5 ' 6 , and exhibit a quite remarkable anomalous elasticity, with wavevector 
dependent elastic moduli that vanish and a bending rigidity that diverges at long wavelength 7 . 
Excluded volume interactions, present in physical membranes, further stabilize the flat phase 8 
but are usually assumed to be otherwise irrelevant to describe its long distance properties. 
Motivated by recent experiments on partially polymerized vesicles 3 , studies of models with 
quenched in-plane disorder have shown that the flat phase is unstable at T=0 to either local 
random stresses 9 or random spontaneous curvature 10 . 

Flat membranes of internal dimensionality D and linear size L are characterized by a rough- 
ness exponent ( such that transverse displacements scale as IS. Nelson and Peliti (NP), using 
a simple one loop self-consistent theory 5 for D = 2 which assumes non-vanishing elastic con- 
stants, found that phonon-mediated interactions between capillary waves lead to a renormalized 
bending rigidity Kn(q) ~ q~ v with rj = 1. Since £ = (4 — D — i])/2 they predicted ( = 1/2 for 
physical membranes. An e = 4 — D expansion 7 confirmed that the flat phase was described by 
a non trivial fixed point, but with anomalous elastic constants \{q) ~ /i(q) ~ q Vu , i] u > 0, with 
i] u = 4 — D — 2i] as a consequence of rotational invariance. Thus, in general ( = (4 — D + r] u )/4 
and the NP approximation corresponds to setting i] u = 0. 

There is presently some uncertainty on the precise value of the roughness exponent for 
physical membranes. Numerical simulations of tethered surfaces display a range of values for 
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( from 0.5 2 , 0.53 11 , 0.64 8 ' 12 to 0.70 13 . On the other hand, the 0(e) result 7 suggests a value 
very close to the NP value 1/2, ( 0.52 by naively setting e = 2). ( should soon be measured 
from experiments, either directly from light scattering on diluted solutions 4 or indirectly, from 
the scale dependence of the elasticity 14 of lamellar stacks of solid membranes presently under 
experimental study. The buckling transition 6 , if observed, is controlled by a single exponent 
related to (. It thus seems desirable to explore further possible theoretical predictions for (. 

In this Letter we introduce a self-consistent approximation which improves on the Nelson 
Peliti theory 5 by allowing a non trivial renormalization of the elastic moduli. It is exact in three 
different limits and compares well with numerical simulations. We construct two coupled self- 
consistent equations for the renormalized bending rigidity Kn(q) and elastic moduli fin(q), \n(q) 
and solve them in the long wavelength limit. Kn(q) is determined by the propagator for the d c = 
d — D components h of the out of plane fluctuations G(q) ~ l/g 4_D ~ ,? while the elastic moduli 
are determined by the four-point correlation function of h fields. Physically, our calculation 
includes the additional effect of relaxation of in-plane stresses by out of plane displacements. 
As a result, curvature fluctuations soften elastic constants and screen the phonon- mediated 
interaction. A similar Self- Consistent Screening Approximation (SCSA) was introduced by 
Bray 15 to estimate the r\ exponent of the critical 0(n) model (here d c plays the role of the 
number of components n) and amounts to a partial resummation of the l/d c expansion. By 
construction, the method is exact for large codimension d c to first order in l/d c and arbitrary 
D. Solving self-consistently then leads to an improved approximation of fj(d c , D) (and thus () 
for the small (physical) values of d c . 

The attractive feature of our theory is that it becomes exact in several other limits. Firstly, 
because of the Ward identities associated to rotational invariance we find that rj(d c , D) is exact 
to first order in e = 4 — D for arbitrary d c and is thus compatible with all presently known 
results 6 ' 7 . Secondly, for d c — it gives rj = (4 — D)/2 which is the exact result since clearly 
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r\ u = for d = D, and 7 r\ u = 4 — D — 2rj. This is at variance with the 0(n) model for which 
the SCSA 15 is not exact for n — 0. Thus we expect this method to give more accurate results 
for the present problem. Two-loop calculations are in progress 16 to estimate the deviation. 
An encouraging indication is the similarity of our method with the remarkably accurate self- 
consistent approximation of Kawasaki 17 for the critical dynamics of the binary fluid mixture, 
which was shown to be exact to order e, again because of Ward identities, and incorrect to 
order e 2 by a tiny amount. We also apply this method to the crumpling transition of phantom 
membranes, and to flat membranes with quenched disorder. Details can be found in Ref.16. 

In the flat phase, the membrane in-plane and out-of plane displacements are parametrized 
respectively by a D-component phonon field u a (x), a — 1, ..D, and a d c = d — D component 
out-of plane height fluctuations field h(x). A monomer of internal coordinate x is at position 
r(x) = (x a + u a (x))e a + h(x) where the e a are a set of D orthonormal vectors. The effective 
free energy is the sum of a bending energy and an in plane elastic energy (most relevant terms): 

F = /d^(W + ^ + ^«L] (i) 

where the strain tensor is u a p = \(d a up + dpu a + d a h ■ dph) To discuss the SCSA in the flat 
phase it is convenient to first integrate out the phonons 1,5 , and to work with the d c -component 
h field. In terms of Fourier components the free energy takes the form of a critical theory: 

F eff = ^ J dkk 4 | h(k) | 2 J dkidk 2 dk 3 R al3 ,js{q) k la k 2 pk 3l k4s h(fci).h(A;2) h(fc 3 ).h(A;4) 

(2) 

with q = k\ + k 2 and k\ + k 2 + k 3 + k± = and we use / dk to denote / d D k/{2n) D . The 
four-point coupling fourth-order tensor R(q) is transverse to q, the longitudinal part having 
been eliminated through phonon integration. It can be written as R(q) = bN(q) + /iM(q) with: 

Nap,* = p^jP&PZ , M aPn& = l -{Pl^ s + Pl s Pl)-N a ^ s (3) 

where P^ = 5 a p — q a q$lq 2 is the transverse projector, /i is the shear modulus and b = 
/i(2/i + DX)/(2/i + A) is proportional to both shear and bulk moduli. The convenience of this 
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decomposition is that M and N are mutually orthogonal projectors under tensor multiplication 
(e.g M^sM^nu = M aP ^ u etc.). 

We set up two coupled integral equations for the propagator of the h field and for the 
renormalized four point interaction. We want to evaluate < hi(—k)hj(k) >= 5ijG{k) with 
G^ik) = n R {k)k A = nk A + a(k) where a(k) is the self energy. The SCSA is defined in 
diagrammatic form by the graphs of Fig. la and lb, where the double solid line denotes the 
dressed propagator G(q), the dotted line the bare interaction R(q) and the wiggly line the 
"screened" interaction R(q) dressed by the vacuum polarization bubbles. We thus obtain two 
equations, one for a{k) which determines 77, the other for R which determines i] u : 

2 f 

cr(k) = —kakpk^ks J dqR aPa& (q)G(k - q) (4a) 

R(q) = R(q) - R(q)U(q)R(q) (46) 

where H a/ 3^s(q) = / dpPaP/3PyP5G(p)G(q — p) is the vacuum polarization and tensor multipli- 
cation is defined above. Because of the transverse projectors, only the component U(q) sym of 
n(g) proportional to the fully symmetric tensor S a p tl s = ^a^-yS + ba-fips + ^a5<Vr contributes in 
(4b). Defining H(q) sym = I(q)S, simple algebra gives R(q) = p,(q)M + b(q)N with renormalized 
shear and shear-bulk moduli, and the new equations: 

m = T^m ~ b(q) = TTini TjTWb (5a) 

We now solve these equations in the long- wavelength limit. Substituting G(k) ~ o" _1 (^) ~ 
Z/k A ~ v in (5a,b), with Z a non-universal amplitude, we find that the vacuum polarization 
integral diverges as: 

I(q) ~ Z 2 A(D,7])q-^ (6) 



5 



where r\ u = 4 — D — 2i] is the anomalous exponent of phonons. Substituting in (5a,b), and 
defining the amplitude: / dqq^(k - q)- ( - 4 -^(kP T (q)k) 2 = B(D, 77)^, one finds (for fi, b > 0) 
that the Z and k 4 ~ v factors cancel and that rj is determined self-consistently by the equation 
for the amplitude: d c = 7^4 a(£>'^] ' which after calculation of the integrals defining A,B gives: 
. _ l D(D _ n T[l + ^]T[2 - V ]T[ V + D]T[2 - \ V ] 

For D = 2 this equation can be simplified, and one finds (Fig. 2): 

V (D = 2, dc) = . 4 (8) 

d c + y/16 - 2d c + dl 

Thus for physical membranes we obtain: rj = 0.821, r/ u = 0.358 and: 

roughly at midvalue of the present numerical simulations. From (5) we also obtain Lim X(q) / jl(q) 
(i.e a negative Poisson ratio). 
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Expanding the result (7) in l/d c one obtains: 

"=^§Ti r [ fi3 r rg- f i + °^ ) 4 + 0( 4 ) (! ° rD = 2) (10) 

which coincides with the exact result 6,7 , as expected by construction of the SCSA. Similarly, 
expanding (7) to first order in e = 4 — D one finds: 

V= 2TdJ^ 

also in agreement with the exact result 6 ' 7 . This is not a general property of SCSA. Here it 
can be traced to the vertex and box diagrams of Fig. lc being convergent. Indeed, because 
of the transverse projectors in (2-3) one can always extract one power of external momentum 
from each external h legs, which lowers the degree of divergence from naive power counting. 
As a result, if one decouples the 4-point vertex R via a mediating field, the only counterterms 
needed are for two-point functions. 
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We have analyzed the crumpling transition of phantom membranes by the same method, 
applied to the isotropic theory of Ref.18. The exponent i] = i] cr at the transition is determined 
by 16 : 

D(D + 1)(D - 4 + ??)(£- 4 + 2r])(2D - 3 + 2^)r[^]r[2 - r]]T[r] + D]T[2 - \r,\ 
~ 2(2 - 77)(5 - D - 2 V )(D + 77 - 1)T[±D + ^]r[2 - 77 - \D]T[r] + \D\V[\D + 2-\rj{ 

(12) 

At the transition the radius of gyration scales as Rq ~ L v with v — (4 — D — rj cr )/2. For 
d = 3 and D = 2 we find ?7 cr = 0.535 and z/ = 0.732 (Haussdorf dimension d H = 2.73). The 
embedding dimension d u (D) above which self-avoidance is irrelevant for the membrane at the 
crumpling transition is determined by the condition d u = 4D /(4 — D — rj cr (d u )). Using (12) we 
find that d u {2) =4.98. 

The present method gives interesting predictions for lower critical dimensions. In the flat 
phase, orientational order (i.e in Vh ) disappears for D < D lc , where 2 — r)(Di c , d c ) = D Xc . 
From (7) this is equivalent to d c = Di c (D lc — l)/(2 — D lc ). On the other hand, the lower critical 
dimension D' lc {d) for the crumpling transition is defined by 2 — r) cr (D' lc , d) = D' lc , or equivalently 
from (12), d = D[J(2-D[ C ). Since d = D/(2-D) is clearly equivalent to d c = D(D-l)/(2-D) 
we find that the lower critical dimensions of the crumpling transition and of the flat phase, as 
predicted by SCSA, are identical, and given by Di c (d) = 2d/{l + d). Since they originate from 
very different calculations, this indicates that the SCSA is quite consistent. For codimension 
1 manifolds Ac = and for fixed embedding space d = 3, Di c = 3/2. Di c increases from 
Di c — 1 for d c — 0, to Die = 2 when d c — > 00 as expected. Note that for d > 3 self- avoidance 
cannot modify the above results, while for d < 3 it is an open question. 

We can compare (8,12) with recent simulations 19 of D = 2 membranes with self-avoidance 
in higher d c . The membranes are found flat in d = 3, 4 with ((d = 3) = 0.64 ± 0.04, ((d = 4) = 
0.77 ± 0.04, whereas we obtain 0.59, 0.67, respectively. The membrane is crumpled in d — 5 



7 



with v = 0.8 ± 0.06, although d — 5 seems almost marginal, whereas we find v = 0.8 at the 
crumpling transition where self avoidance is irrelevant, although almost marginally so. 

Flat membranes with random spontaneous curvature are described by adding the term 
— / d D 'xc(x) .V 2 h(x) in the energy (1), where c(x) are Gaussian quenched random variables 10 . 
Within a replica symmetric SCSA, we find a marginally unstable T = fixed point, i.e a 
long- wavelength solution only if T — > first. Defining the replica connected and off-diagonal 
exponents 77, r/, by < h(— q)h(q) > c ~ q~( 4 ~ v \ < h(— q)h(q) > ~ q-i^-v') we find 16 at this 
fixed point: rj' = rj, rj(d c , D) = i] pure (4:d c , D). Thus one can simply replace d c in the pure result 
by Ad c \ Again this agrees with the l/d c and e expansions 10 . For physical membranes D = 2, 
d c — 1, we find from (8) : 

r] = 2/(2 + VQ) = 0.449 C = 0.775 

comparing well with the numerical simulation 10 result ( = 0.81 ± 0.03. By analogy with the 
random field problem 20 , it is quite possible that the equality i] = rf ', conjectured in Ref. 10 to 
all orders, be corrected when replica-symmetry breaking is included. 

In conclusion, we have presented a self-consistent theory of polymerized membranes which 
becomes exact in three limits (large d c , small e = 4 — D, and d c = 0). By construction, it 
satisfies the exponent relations rj u = 4 — D — 2i] and 16 — D — 2 + 77. These relations are 
exact in the true theory because of rotational invariance 6 ' 7 . It thus predicts v' = 1.218 and 
5' = 1.436 for the buckling transition exponents 6 . It contradicts the conjecture 2 ( = 1/2. 

We thank D. Nelson, M. Mezard for discussions. PLD acknowledge support from NSF grant 
DMS-9100383 and LR from the Hertz Graduate Fellowship. 



8 



References 

*Also LPTENS, Ecole Normale Superieure, 24 rue Lhomond, Paris 75231 Cedex 05, Labo- 
ratoire Propre du CNRS. 

1. See, e.g, Statistical Mechanics of Membranes and Interfaces, edited by D.R. Nelson, T. 
Piran, S. Weinberg ( World Scientific, Singapore 1988 ), and S. Leibler in Proceedings of 
the Cargese school on biologically inspired physics, (1990), to be published. 

2. R. Lipowsky and M. Girardet, Phys. Rev. Lett. 65 2893 (1990). 

3. M. Mutz, D. Bensimon, M.J. Brienne, Phys. Rev. Lett. 67 923 (1991). 

4. X. Wen et al. Nature 355, 426 (1992) 

5. D. R. Nelson and L. Peliti, J. Phys. (Paris) 48, 1085 (1987). 

6. F. David and E. Guitter, Europhys. Lett. 5, 709 (1988). E. Guitter, F. David, S. Leibler 
and L. Peliti, J. Phys. France 50 1787 (1989). 

7. J.A. Aronovitz and T.C. Lubensky, Phys. Rev. Lett. 60, 2634 (1988), J.A. Aronovitz, 
L. Golubovic, T.C. Lubensky, J. Phys. France 50 609 (1989). 

8. F.F. Abraham, D.R. Nelson, J. Phys. France 51 2653 (1990). F.F. Abraham, W.E. Rudge 
and M. Plishke, Phys. Rev. Lett. 62, 1757 (1989). 

9. D.R. Nelson, L. Radzihovsky, Europhys. Lett. 16, 79 (1991), L. Radzihovsky, P. Lc 
Doussal, J. Phys. I France 2, 599 (1992). 

10. D.C. Morse, T.C. Lubensky, G.S. Grest, Phys. Rev A 45 R2151 (1992). Morse Lubensky, 
Preprint 1992. 



9 



11. F. Abrahams, Phys. Rev. Lett. 67, 1669 (1991). 

12. S. Leibler and A. Maggs, Phys. Rev. Lett. 63, 406 (1989). 

13. G. Gompper and D.M. Kroll, J. Phys. I France 2, 663 (1992). (1992). 

14. J. Toner, Phys. Rev. Lett. 64 1741 (1990). 

15. A.J. Bray, Phys. Rev. Lett. 32, 1413 (1974). 

16. P. Le Doussal, L. Radzihovsky, to be published. 

17. E.D. Siggia, B.I. Halperin and PC. Hohenberg, Phys. Rev. B 13 2110 (1976). 

18. M. Paczuski, M. Kardar and D.R. Nelson, Phys. Rev. Lett. 60, 2638 (1988). 

19. G. Grest, J. Phys. I France 1,1695 (1991). 

20. M. Mezard, A.P Young, Preprint LPTENS 92/2 



10 



FIGURE CAPTIONS 



Figure 1 : graphical representation of the SCSA: (a) self energy, (b) interaction, (c) UV finite 
vertex and box diagrams. 

Figure 2 : ( as a function of d for two-dimensional membranes D = 2. Solid curve: SCSA 
result (8). Dashed dotted curve: 0(e) result, setting e = 2. Dashed curve: corresponds 
to 7] = 2/d chosen (somewhat arbitrarily) in Ref. 6 as a possible interpolation to finite d 
(asymptotic to the solid curve for d — > oo). 
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